Scattering Theory for Quantum Hall Anyons in a Saddle Point Potential 
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We study the theory of scattering of two anyons in the presence of a quadratic saddle-point 
potential and a perpendicular magnetic field. The scattering problem decouples in the centre-of- 
mass and the relative coordinates. The scattering theory for the relative coordinate encodes the 
effects of anyon statistics in the two-particle scattering. This is fully characterized by two energy- 
dependent scattering phase shifts. We develop a method to solve this scattering problem numerically, 
using a generalized lowest Landau level approximation. 
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One of the most remarkable features of the fractional 
quantum Hall effect P, [3] is the existence of quasipar- 
ticles with fractional charge Q and fractional exchange 
statistics 0, 0. Direct evidence for the existence of frac- 
tional charge has been observed experimentally 0,0, 
One way this has been achieved is through the study of 
shot noise in a quantum Hall system constricted by a 
point contact [?', '«|, by which the charge of the current 
carriers may be derived from the fluctuations of the cur- 
rent backscattered from the point contact. 

Despite strong theoretical reasons to expect that these 
fractionally charged quasiparticles also have fractional 
statistics [j, to date there has been no unambigu- 
ous experimental demonstration of fractional statistics. 
Several authors have proposed methods in which evi- 
dence of fractional statistics might appear in transport 
experiments [lol. [ll . 1^ 13, 13], based on mesoscopic de- 
vices of various geometries. Here we adopt a simple ap- 
proach, and ask a very natural theoretical question: can 
anyonic statistics influence the transport through a sin- 
gle point contact? Clearly, since the effect of exchange 
statistics is a two-particle property, any such effect re- 
quires more than one particle to be present in the point 
contact region. We therefore study the scattering prop- 
erties of a pair of anyons incident on the point contact. 
Since the quantum Hall anyons are charged, and experi- 
ence a strong magnetic field, their free motion is along 
edge states which follow classical equipotentials. The re- 
quired scattering theory is therefore very different from 
that of freely moving anyons, described by a conventional 
kinetic energvfisj]. As we show, the anyonic nature does 
affect the scattering of two quantum Hall quasiparticles 
in the point-contact region. This effect can be charac- 
terized by two energy-dependent phase shifts, which we 
calculate numerically for general anyonic statistics pa- 
rameter. 

We consider two identical particles, with mass M and 
charge q, subjected to a magnetic field B and in the pres- 
ence of a quadratic saddle point potential 



H 



^ 2M 

i=l,2 



where V x A = Bz is the uniform magnetic field. The 
particles are taken to be anyons [l6|, with statistics pa- 
rameter A. Thus, the two-particle wave function has the 
boundary condition that 



27r) = e 



-i27rA 



(2) 



where 9 ^ 9 + 2tt represents one complete clockwise ro- 
tation of the relative co-ordinate of the two particles. We 
note that this defining feature applies either for distin- 
guishable or indistinguishable anyons 17J. If the anyons 
are indistinguishable, one can impose the stronger re- 
striction 

^' (61 + tt) = e-^'''^-^ (9) [indistinguishable] (3) 

Then, one recovers the exchange statistics for bosons and 
fermions for A = 0, 1 respectively. For greater generality, 
we study the case of distinguishable anyons. The effects 
of the additional restriction for indistinguishability ([3]) 
will be made clear in the discussion. 

We simplify the problem, defined by ^ and the 
boundary condition (l2|), by introducing a "centre-of- 
mass" coordinate, Vc, and a "relative" coordinate, r^, 
with 



= (ri -I- rz) Tr = (ri - rs) 
The Hamiltonian separates, becoming 



Ha 



fj,—c,r 



(4) 



(5) 



qA{r,)f + U {yf ~ xj) (1) 



Thus, the two-particle problem can be expressed as two 
independent one-particle problems, and the total energy 
E = Ec + Er can be divided into separately conserved 
contributions from the centre-of-mass and relative co- 
ordinates. 

The centre-of-mass co-ordinate is insensitive to any- 
onic statistics. The scattering theory for this co-ordinate 
is identical to the one-particle problem solved by Fertig 
and Halperin[3]. This work allows one to deduce the 
transmission and refiection coefficients of the centre-of- 
mass co-ordinate moving in the saddle point potential, in 
terms of Er. 
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The relative coordinate has the same Hamiltonian. 
However, since the wave function has the additional any- 
onic boundary condition the results of Ref.[l8|] can- 
not be applied. The solution of the scattering problem 
for the relative co-ordinate is the central result of the 
present paper. From here on, for simplicity, we drop the 
subscripts r on and Er, with the understanding that 
the calculation refers only to the relative co-ordinate. 

Before considering the scattering problem, we study 
first the spectrum for the relative co-ordinate in the ab- 
sence of a potential, U = 0. For A = 0, this leads to the 
familiar Landau level states and spectrum. The anyon 
boundary condition changes the nature of these states. 
To describe the generalized Landau level states, we use 
polar coordinates, r = r{sin9,cos9) and the symmetric 
gauge A = —eorB/2. The eigenstates are 



_ g-i(m+A)6 



(6) 



where m is an integer, and n — 0, 1, 2, . 
ized radial wave functions are 



The normal- 



(7) 



where £ = y^h/ {qB) is the magnetic length, are 
the Laguerre polynomials, and the normalization is Af = 

^nl/ {r{\m + A\ + 



1) 27r2l™+^l). The energy is 



1 



-|m + A|-i(m + A) + i 



(8) 



where u>c = qB/M is the cyclotron frequency. The spec- 
trum has the required feature that it is invariant under 
the transformation A ^ A' = A — 1, equivalent to the 
insertion of a flux quantum at the origin. This amounts 
to the change m — s- m' = to -I- 1 . Thus, it is sufficient to 
study the range < A < 1 to cover all possible cases. 
For positive to -|- A > the energy depends only on n: 
these sets of states thus form highly degenerate Landau 
levels for the relative motion of the anyons. 

We now reintroduce the potential U ^ 0. We do this 
within a (generalized) lowest Landau level approxima- 
tion, in which U is taken only to lift the degeneracy 
of the lowest energy states. This corresponds to re- 
taining only those levels with n — Q and to -1- A > 0. 
Noting that the lowest energy state lies higher in en- 
ergy by ?iu;c(l — A), the lowest Landau level approxi- 
mation for the relative motion of the anyons is valid for 
^ fiLOc{l ~ A). We denote the set of degenerate 
basis states as |j) where j = 0, 1, 2, . . ., and expand the 
wave function as jV') = X^V'ib)- The Schrodinger equa- 
tion becomes 



^{j + A - 1) (j + A)V;,_2 
+ v/(j + A + 1) (j + A + 2)^, 



J+2 



(9) 



for the amplitudes V'j, where e = E/iUi"^) is the dimen- 
sionless measure of the energy. Note that there are no 
terms that couple odd and even values of j. Physically, 
this arises from the fact that the potential is invariant un- 
der spatial inversion {x, y) — > {—x, —y), and so the parity 
of the wave function is a good quantum number. Thus 
the Schrodinger equation takes the form of two decou- 
pled sets of difference equations. We write the general 
solution in terms of the "even" and "odd" channels as 



p=0 



p=0 



(10) 
(11) 



For distinguishable anyons, the wave function can be a 
linear superposition of these two solutions. However, for 
indistinguishable anyons, the boundary condition re- 
quires that only the even solution contributes. 

We shall construct the wave function at large distances 
from the origin, r ^ i. In this limit, contributions are 
from single particle states with j ^ 1. For j ^ e, 1, 
the Schrodinger equation ^ has the wave-like solutions 
(normalized to unit density per orbital j) [l9| 



[00 = ±7r/4,±37r/4] 



(12) 



These solutions can be viewed as incoming and outgoing 
waves in the discrete semi-infinite one-dimensional sys- 
tem defined by the sites j — 0, 1, 2 . . .. To understand 
the nature of these solutions, it is useful to construct 
their spatial wave functions 



(r|V) = J2 



"^^■+^^'i?o,, (r) 



(13) 



At large radius r, the wave function has significant am- 
plitude under the condition that ~ 0o = i7''/4, ±37r/4. 
That is, the wave function is peaked in these four angu- 
lar directions. These are the directions along which the 
zero energy equipotentials of the electrostatic potential 
extend. Recalling that, in the semi-classical approxima- 
tion, the particle moves along the equipotentials of the 
electrostatic potential, one sees that the derived angles 
correspond to the two incoming (— 7r/4, 37r/4) and the 
two outgoing (7r/4, — 37r/4) channels of the saddle-point 
potential. 

For a general scattering problem on a semi-infinite 
one-dimensional system (e.g. waves on a string which 
is clamped at one end) , at a fixed energy (frequency) one 
expects there to be only two wave-like solutions at large 
distances; these can be taken to be the incoming and the 
outgoing waves. That, in the present case, there are jour 
wave-like solutions is a special feature of the problem, 
which arises from the fact that (as above) the Schrodinger 
equation ([5]) conserves the parity. That is, the sites with 
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j even and the sites with j odd each behave as inde- 
pendent semi-infinite one-dimensional systems. For each 
parity (even or odd) there is one incoming mode and one 
outgoing mode. It is convenient to re-express the four 
modes ([T^ in terms of modes of definite parity. The 
exphcit forms (normahsed to unit density per orbital j) 



Noting that the transmitted wave is 



are 



c/o,in 



c/o,out 



1 

V2 

1 

V2 



g»f j _|_ g-«TJ 



(14) 
(15) 



which are readily verified to have the feature that -0° 
is non-zero only for j = even (odd). The "in" and "out" 
labels are identified by the fact that the states have large 
amplitude on the incoming (6 — — 7r/4, 37r/4) or outgoing 
{6 = 7r/4, — 37r/4) channels of the saddle point potential. 

For the semi-infinite one-dimensional scattering prob- 
lem, the asymptotic (large distance) incoming and out- 
going waves are coupled, with scattering from incoming 
to outgoing waves occurring at small distances. In the 
problem of interest here this scattering at small distances 
(small j) conserves parity, so the incoming mode in the 
even (odd) channel can scatter only to the outgoing mode 
in the even (odd) channel. By conservation of particle 
fiux, the scattering from incoming to outgoing modes can 
amount only to a phase shift. Hence, the (mmormalised) 
energy eigenstates are of the form 



e/o 



, c/o.in 



3''C^°(e)^<=/o,out 



(16) 



Thus, the energy eigenstates are fully characterized by 
two scattering phases: C° (e) and (° (e) which correspond 
to the even- and odd-parity wave functions. 

For indistinguishable anyons ([3]) only is relevant. 
There is only one incoming and one outgoing channel 
for the relative co-ordinate, so the scattering is described 
only by a single phase shift, C°(e)- 

For distinguishable particles, both -0° and ip° can con- 
tribute. In this case, it is instructive to disentangle the 
above transformation into states of definite parity, and 
to determine the transmission probability. This is defined 
as the probability for transmission from a state that is 
an incoming wave along the definite angular direction 
6 ~ — 7r/4, into a state that is outgoing along 9 ~ -|-7r/4. 
(These angles match the convention chosen in Ref. [l8|.) 
The incoming wave (normalised to unit density per or- 
bital) is 



U - 



1 

V2 



which is scattered into the state 



1 



/ c.out , iC" I o.out 

V'/ +e V/ 



(17) 



(18) 



4--i4 7 1 r / c.out 



(19) 



and using the fact that and are orthogonal, 

one sees that the transmission probability is 



T 



+ e 



1 



{l + cos[C {e)-C (e)]} 



(20) 

We have determined the functions C'^^°{e) - thereby 
solving the scattering problem for the relative co-ordinate 
- by a numerical construction of the Green's function of 
the discrete Hamiltonian in ([9]). We study an approxi- 
mate version of the full model, in which we treat j = 0, N 
according to the exact Hamiltonian ([9]) , but take the hop- 
ping matrix elements for j = + 1, cx3 to be constant, 
and equal to ^/{N + A - 1) {N + A). The method be- 
comes increasin gly accurate as ^ cx). Following stan- 
dard techniques [205, the region with j > -I- 1 can be 
replaced by a self-energy, and the Green's function for 
j < N is: 



(21) 



G{e) = el- H -a 



where the self-energy has matrix representation 

^■u = I^nSj^n (e - iV4(iV + A+ 1) (^ ~t- A + 2) - 

(22) 

The Green's function for j < N is then found by numer- 
ical inversion of the finite matrix (for large finite N). 

Using the fact that, for a given element i, the Green's 
function Gij for j < i is an energy eigenfunction, we 
can use this to compare to the wave function (|16[) in the 
asymptotic regime j ^ 00. Using (fT6|) . the scattering 
phase can be extracted from the ratios 



■ip+2 P- 



i'lp+l 



tan 



tan 



(e) + I ln(4p) 



-C° (e) - 7r/2 -f f ln(4p) 



(23) 



(24) 



where the logarithms on the right-hand side follow from 
the corrections to the wave functions (fT^ described in 
[l9| . In this way we can numerically construct the 
scattering phases C°/° (e) and from these the transmis- 
sion probability T (|20p . The numerical results converge 
rapidly with increasing iV, becoming independent of sys- 
tem size for N > 100. We show results for N = 2000. 

Figs. [1] and [2] show our numerical results for the scat- 
tering phase shifts for the even and odd parity channels 
respectively, as a function of energy for several values 
of A. These two functions, over the range < A < 1, 
fully describe the scattering properties of the relative co- 
ordinate of quantum Hall anyons in the lowest Landau 
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FIG. 1: (Colour online.) Results of numerical calculations 
of the scattering phase shift for the relative motion of two 
anyons in the even parity channel, C^i^), as a function of di- 
mensionless energy e — E/Ul^ at several values of the anyonic 
statistics parameter A. 
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FIG. 2; (Colour online.) Same as Fig. [T] but for the odd 
parity channel, C°(^)- 



level. The results shown for A = 1 correspond to the 
case A = 1— in which the state m = — 1 remains ex- 
cluded from contribution to the lowest Landau level. In 
this case, the spectrum is identical to that for A' = 
and m! = m -f 1, but with the removal of the state at 
to' = 0. From (I16p . and taking account of a 7r/2 phase 
shift arising from the change m' = to + 1, one expects 
CA=i_(e) = Ca'=o(^) + ""/^j which is indeed found to 
hold to high accuracy in the numerical results. In Fig. [3] 
we have used these results to determine the transmission 
coefficient (|20]) . For A = 0, the results accurately re- 
produce the exact analytical solution[l^, showing that 
our method is working correctly and is well converged. 
For A 7^ 0, the results of Ref.[l8| do not apply, and no 
analytic solution is available. As compared to the case 
A = 0, the effect of increasing A is a broadening of the 
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FIG. 3: (Colour online.) Transmission probability T for the 
relative coordinate of two quantum Hall anyons, as a function 
of the scaled energy e = Er/{Ul'^). The points are results of 
numerical calculations for several values of the anyon statistics 
parameter A. The lines are fits, using the function (|25|l with 
parameter a chosen as shown in the inset. For A = choosing 
OL = 1 gives the exact analytic result in this case, the 
comparison with numerical results illustrates the accuracy of 
our numerical approach at this system size (A'' = 2000). 



width in energy over which the transmission coefficient 
rises from (e^ —1) to 1 (e^ 1). Thus, our results 
show that increasing A leads to an increase of the tun- 
nelling rate through the saddle point. Indeed, we find 
that, to a very good approximation, the results can be 
fitted by the function 



r(e) = 



1 



1 



exp I 



-) 

OL } 



(25) 



with the fit parameter a given approximately by a = 
1 + 1.55378A ^ 0.277179A2. (See the inset to Fig. [3) 

We note that our results apply for the case of two 
anyons in a symmetrical saddle point potential ([T]). The 
more general case can be considered by noting that 
[7,2/2 _ [/^a;2 = 1 ([/^ + C/,) _ + \{Uy - C/^)r2, 

leading to an additional central (rotationally invariant) 
term cx r^. This term modifies the Schrodinger equa- 
tion (lU and could lead to a change in the phase shifts 
and transmission probabilities. In the same way, central 
(rotationally invariant) anyon-anyon interactions could 
be included within the same formalism. The solution 
of these more general cases is beyond the scope of the 
present paper, so the influence of these perturbations on 
the scattering properties remains an open question. 

In summary, we have provided a solution of the scat- 
tering problem for two anyons in a quadratic saddle point 
potential and perpendicular magnetic field, through sepa- 
ration into centre-of-mass and relative coordinates. The 
scattering theory for the centre-of-mass coordinate has 
previously been solved analytically [3]. The scatter- 
ing for the relative co-ordinate is characterized by two 
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energy-dependent phase shifts (for even and odd par- 
ities). We have computed these phase shifts within a 
lowest Landau level approximation. Our results provide 
a complete solution of the two-anyon scattering prob- 
lem. They show that the two-particle scattering prop- 
erties in the vicinity of a point contact depend on the 
anyonic statistics parameter. This shows that one can 
hope to obtain experimental signatures of anyonic statis- 
tics in point-contact devices. We believe that the ap- 
proach and solution we have described provide a useful 
basis on which to build further theoretical studies of non- 
equilibrium properties of anyons - under conditions of 
bias and/or "beam dilution" [2lj where multiple quasi- 
particles may enter the point contact region. Depending 
on the experimental set-up, and the relevant observables, 
the two-particle scattering problem in the co-ordinates 
ri , r2 should be decomposed into relative and centre-of- 
mass co-ordinates for which the scattering properties fol- 
low the results presented. 

NRC acknowledges the support of EPSRC 
GR/S61263/01. 
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